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Liquid-state NMR quantum computer has demonstrated the possibility of quantum computation 
and supported its development. Using NMR quantum computer techniques, we observed phase 
decoherence under two kinds of artificial noise fields; one a noise with a long period, and the other 
with shorter random period. The first one models decoherence in a quantum channel while the 
second one models transverse relaxation. We demonstrated that the bang-bang control suppresses 
decoherence in both cases. 
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I. INTRODUCTION 

Quantum computation currently attracts a lot of at- 
tention since it is expected to solve some of computa- 
tionally hard problems for a conventional digital com- 
puter |l| . Numerous realizations of a quantum computer 
have been proposed to date. Among others, a liquid-state 
NMR (nuclear magnetic resonance) quantum computer 
is regarded as most successful. Demonstration of Shor's 
factorization algorithm J2| is one of the most remarkable 
achievements. 

Although the current liquid-state NMR quantum com- 
puter is suspected not to be a true quantum computer 
because of its poor spin polarization at room temperature 
P|, it still works as a test bench of a working quantum 
computer. Following this concept, we have demonstrated 
experimentally using an NMR quantum coinputer that 
some of theoretical proposals really work P, |^ . In this 
contribution, we will show that a liquid-state NMR quan- 
tum computer can model not only a quantum computer 
but also the composite system of a quantum computer 
and its environment. Therefore, one can employ it to test 
the effectiveness of proposed decoherence control meth- 
ods, such as a bang-bang control 6, 7]. Note that the 
decoherence control methods are usually difficult to be 
tested because of extremely short coherence time in the 
real system. 



Hamiltonian 



Ht — Hs + He 



Hf,, 



(1) 



where Hg and He, in the absence of Hse, determine the 
system and the environment behavior, respectively. On 
the other hand, Hsc determines the interaction between 
the system and the environment. See, Fig. ^ Zurek dis- 
cussed a simplified model where a two-level system (the 
system) is coupled to n two-level systems (the environ- 
ment) through az ^ a^ type interaction Q . 



A. Artificial Decoherence 

If the effect of H^c is small enough to be ignored com- 
pared with those of Hs and H^ in a certain time scale r, 
H^ can be considered as 



TJg — H\ + H2 + H\ 



(2) 



H\ and H^, if Hyi does not exist, determine the behav- 
ior of subsystem 1 and 2, respectively, while H\2 de- 
termines the interaction between the subsystems. See, 
Fig. n Therefore, we may regard the subsystem 1 (2) as 
the system (environment) in the time scale r and that 
the dynamics of the subsystem 1 can model that of a 
certain system. Zhang et al. experimentally studied the 



II. DECOHERENCE 

Decoherence is a phenomenon in which a quantum sys- 
tem undergoes irreversible change through its interaction 
with the environment. This is analyzed using the total 
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FIG. 1: System and environment. The system consists of 
subsystems 1 and 2. 



behavior of ^'^C-labeled trichloroethane, which has three 
spins, using NMR techniques Q and claimed that they 
studied the decoherence. 

We, however, note that a large number of degrees of 
freedom of the subsystem 2 is necessary to observe a 
decoherence-like behavior in the subsystem 1. This con- 
dition is not satisfied with molecules employed in liquid- 
state NMR quantum computation. If the degrees of free- 
dom of the subsystem 2 is small, a periodic behavior in 
the subsystem 1 should be observed instead of an irre- 
versible one. Teklemariam et al. introduced a stochastic 
classical field which is acting on the subsystem 2 in or- 
der to overcome the limitation of the model caused by the 
small degrees of the freedom of the subsystem 2 [l^ . The 
approach by Teklemariam et al. can be considered as a 
generation of random noise on the subsystem 1 through 
the subsystem 2 by applying a stochastic classical field 
to the subsystem 2. 



B. Artificial Decoherence in One-Qubit 

Let us consider a molecule containing two spins 
(qubits) as a system. The first qubit is regarded as the 
subsystem 1 and the second qubit as the subsystem 2. 
The Hamiltonian, when an individual rotating frame is 
assigned to each qubit, is 



H. 



Jh 



(3) 



where Ik = crfe/2 and au is the fc-th Pauli matrix. Note 
that Hi = H2 = in this rotating frame. We take a se- 
ries of TT-pulses acting on the second qubit as a classical 
field introduced by Teklemariam et al. |l3|- We cre- 
ate a pseudo-pure state |00) before starting experiments. 
Therefore, the Hamiltonian ^ is equivalent with 



H, 



J{t)h. 



(4) 



The spin operator I^ appeared in Eq. |0J denotes the spin 
of the subsystem 1, while \J{t)\ — J and its sign changes 
when the 7r-pulse acts on the spin 2. We assume that 
the duration of a 7r-pulses is infinitely short. Therefore, 
we can model "a system containing one qubit in a time 
dependent field" , where the field strength is constant but 
its sign changes in time (telegraphic). 

A stochastic classical field is, here, a series of 7r-pulses 
acting on the second qubit randomly in time. 



III. EXPERIMENTAL SET-UP 

A 0.6 ml, 200 mM sample of ^^C-labeled chloroform 
(Cambridge Isotope) in d-6 acetone is employed as a two- 
qubit molecule and data is taken at room temperature 
with a JEOL ECA-500 NMR spectrometer, whose hy- 
drogen Larmor frequency is approximately 500 MHz |llj | . 
The measured spin-spin coupling constant is J/27r — 
215.5 Hz and the transverse relaxation time is T2 ^ 7.5 s 



for the hydrogen nucleus (subsystem 2) and T2 ~ 0.30 s 
for the carbon nucleus (subsystem 1). The longitudinal 
relaxation time is measured to be Ti ^ 20 s for both nu- 
clei. The duration of a 7r-pulses for both nuclei is set to 
50 /^s. 



A. Decoherence in Channel 

Let us consider a flying qubit traveling in a channel, 
as a first example. It is assumed that there exists a noise 
source on a certain position in the channel, which causes 
decoherence in the flying qubit. 

In order to model the above case, we performed an 
experiment schematically shown in Fig.|21 The pseudop- 
ure state loo) (or, |tt)) is prepared by the field gradient 
method [IJ. A 7r/2-pulse acts on the spin 1 aX t — Q. 
Then, the spin 1 is turned into the a;-axis in the rotating 
frame and starts rotating in the xy-plane with an angular 
velocity —J/2. A pair of 7r-pulses are applied to the spin 
2 at t = ti and ti+S and then the spin 1 rotates with the 
angular velocity J/2 during this period. The FID (free 
induction decay) signal at t = tm is measured and shown 
as the open square in Fig. |21(a). The x- (y-) component 
is the real (imaginary) component of the FID signal. The 
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FIG. 2: Model of "decoherence in a channel", (a) Tiie dy- 
namics of the spins is schematically shown. The spin 1 ro- 
tates in the a:y-plane with the angular velocity —J/2 except 
between t — ti and t = ti+S, during which the angular veloc- 
ity is J/2. The spin 2 is flipped at i = ti and t = ti + 5. The 
FID signal of the spin 1 is measured at t — tm- (b) The pulse 
sequences for realizing the spin dynamics shown in (a). The 
short bar indicates a 7r/2- pulse acting on the spin 1, while 
the long squares are vr-pulses acting on the spin 2. (c) Pulse 
sequences compensating the "decoherence" obtained in (b). 
The 16 long bars are vr-pulses acting on the spin 1 with the 
interval of 0.3 ms. 
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FIG. 3: Experimental results of "decoherence in a channel". 
The FID signals measured at t = tm in Fig. |5| are shown 
in the xy- (real and imaginary)-plane. The solid and open 
symbols are with and without the decoherence suppression 
pulses shown in Fig. |5|(c). (a) The open squares denote the 
first 8 FID signals without averaging. The 8 open triangles 
denote the averaged FID signals over 16 experiments. The 
open circle near the origin denotes the averaged value of all 
(16 X 8 = 128) the FID signals. The region marked by the 
rectangle is enlarged in (b). (b) The solid squares denote the 
first 8 FID signals without averaging. The 8 solid triangles 
denote the averaged FID signals over 16 results. The solid 
square on the x-axis shows the averaged value of all (128) the 
FID signals. 



B. Transverse Relaxation 

Let us model a phenomenon called a transverse relax- 
ation next. Suppose that there is a spin in a magnetic 
field (0, 0, Bq). The spin points the z-direction in thermal 
equilibrium. Then, let us turn the spin in the xj/-plane 
by a 7r/2-pulse. The spin starts rotating in the xj/-plane 
with the angular velocity loq = jBq, where 7 is the gyro- 
magnetic ratio of the spin. This spin rotation is called a 
precession and can be observed as a FID signal in NMR. 
If there is no relaxation mechanism, it precesses forever 
and thus the FID signal does not decay in time. The 
transverse relaxation is a phenomena that the spin is still 
in the xy-p\ane but its FID signal decreases in time. The 
transverse relaxation is modeled, in the simplest case, 
as a random walk process 0, UM O"^ ^^^ circle shown in 
Fig.El(a). 

We performed an experiment schematically shown in 
Fig-E The pseudopure state |00) is prepared by the field 
gradient method [i3|- A 7r/2-pulse is applied to the spin 
1 at i = 0. Then, the spin 1 is turned into the a;-axis 
in the rotating frame and starts rotating with an angu- 
lar velocity J. The pulse sequence shown in Fig. 0] (a) 
provides a reference, which is necessary because the in- 
trinsic relaxation cannot be avoided in the experiments. 
The FID signal is measured ai t — tm after a series of 
TT- pulses acting on the spin 2, of which interval is fixed 
to A = 2 ms here. Note that the number of the 7r-pulses 



signal is normalized so that the data point should be on 
the point (1, 0) when 5 = 0. Note that open squares are 
on a circle of unit radius, but with different phases be- 
cause of different duration 5. We randomly choose 128 5s 
between and 2-k/ J and average FID signals, as shown 
in Fig. O (a). Averaging over all signals gives a smaller 
averaged FID signal than that of each FID signal, which 
indicates that decoherence occurred because of the noise 
source in the channel. 



Kitajima, Ban, and Shibata discussed a method to sup- 
press the above decoherence [l^- They argued that a 
series of 7r-pulses acting on the spin 1, while the spin 1 is 
under the influence of the noise source, should suppresses 
the above decoherence. Their idea is essentially the same 
as the field inhomogeneity compensation using the spin 
echo method [lj|- Figure |21(c) shows the pulse sequence 
realizing their decoherence suppression proposal. Since 
we do not know the exact position of the noise source in 
advance, we apply many (16 here) 7r-pulses to the spin 1. 
If a noise source exists within this period (equivalently, 
region when the flying qubit is really moving) of the 16 
TT-pulses, the effect of the noise is greatly suppressed. We 
observed this behavior in our experiments, as shown in 
Fig. O (c). The amplitude of each FID signals are the 
same and the variation of the phases in the ccy-plane is 
remarkably decreased as shown in Fig.|3|(c). Therefore, 
it is clearly seen that decoherence is greatly suppressed. 
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FIG. 4: Pulse sequences for modeling "transverse relaxation" . 
The pseudopure state |00) is prepared first. The spin 1 is 
turned into the a::-axis by the 7r/2-pulse at f = 0. (a) The 
pulse sequence provides a reference. A series of the 7r-pulses 
acts on the spin 2. The interval of the pulses is fixed to 
A = 2 ms. The sequence shown here yields the solid squares 
in Fig. 1^ (b) The pulse sequences to realize the artificial 
transverse relaxation. The intervals between the pulses are 
randomly modulated. See, the text. The artificial transverse 
relaxation is obtained by averaging over various (128) series of 
TT-pulses on the spin 2. (c) The pulse sequence for realizing the 
bang-bang control. In addition to the pulse sequence shown 
in (b), TT-pulses act on the spin 1, of which intervals are 0.5 ms 
< A. 
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FIG. 5: Experimental results demonstrating artificial "trans- 
verse relaxation" and its suppression. The solid squares show 
the amplitude of the FID signals of the spin 1 while the tt- 
pulses are acting on the spin 2 every 2 ms. Those provide a 
reference for the other experimental results. The open squares 
show the amplitude of the averaged FID signals of the spin 
1. We averaged 128 FID signals with various random series 
of TT-pulses acting on the spin 2, here. The open circles show 
the amplitude of the averaged FID signals as in the case of 
the open squares but with the "bang-bang" control. The bro- 
ken (solid) line is the least square fit of a function e~*' ■^^ and 
T2 = 128 (33) ms to the solid (open) squares. 



on the spin 2 determines the period while the spin 1 is 
under the influence of time dependent field, see Eq. Q . 
The amplitudes of the FID signals with various periods 
are shown as the solid squares. When we obtain a faster 
relaxation than this reference, then we can claim that the 
artificial relaxation is realized. 

The artificial transverse relaxation is realized with the 
pulse sequence shown in Fig. 21(b). In contrast to (a), the 
intervals between the pulses are randomly modulated as 
A{l+a Nd), where a is a parameter defining the strength 
of the relaxation and Nd is a variable which obeys the 
normal distribution. We set a = 0.25 in Fig. The 
artificial transverse relaxation is observed when the FID 
signals with various (128 in this experiment) series of tt- 
pulses on the spin 2 are averaged. The open squares in 
Fig. El show the amplitudes of the averaged FID signals 
with various periods. We observe that the relaxation is 
faster than the reference case. 

The pulse sequence to realize the bang-bang control 
[a, I3 is shown in Fig. ^ (c). In addition to the pulse 
sequence shown in Fig.^(b), regular 7r-pulses act on the 
spin 1, of which intervals are 0.5 ms < A = 2 ms. The 
relaxation with the bang-bang control is clearly smaller 
than that without it, as shown in Fig. Therefore, we 
conclude that the effectiveness of the bang-bang control 



is confirmed. 

The experimental results, when the parameter a is 
changed, is shown in Fig. (HI Note that the strength 
of the relaxation can be controlled by changing the 
parameter a. We suspect that there is a non-Markovian 
behavior (non-exponential decay of the amplitude of the 
FID signal) when a = 0.10. We plan to investigate this 
behavior further in near future. 
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FIG. 6: Experimental results demonstrating artificial "trans- 
verse relaxation" with various a values. The open squares, 
open circles, open triangles, and inverted open triangles show 
the results with a = 0.25,0.20,0.15,0.10, respectively. The 
solid squares show the reference as in Fig. |K1 



IV. CONCLUSION 

We generated artificial decoherence (relaxations) using 
liquid-state NMR quantum computer techniques. The 
first type of decoherence takes place in a quantum chan- 
nel, while the other is a transverse relaxation. Then, the 
bang-bang control is applied to the qubit which carries 
a quantum information and we confirmed that it indeed 
suppresses the two types of decoherence. Moreover, we 
have shown that the nature of the decoherence can be 
controlled by changing parameters. 

The artificial decoherence thus generated is still simple, 
but we can extend our approach further. We believe that 
well controlled artificial decoherence will help to under- 
stand various types of decoherence in the real world and 
to develop methods to overcome them towards physical 
realization of a working quantum computer. 

Extended version of this article with detailed theoret- 
ical analysis is in progress and reported elsewhere jig. 
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